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We present an insightful "derivation" of the Langevin equation and the fluctuation dissipation 
theorem in the specific context of a heavier particle moving through an ideal gas of much lighter 
particles. The Newton's Law of motion {mx = F) for the heavy particle reduces to a Langevin 
equation (valid on a coarser time scale) with the assumption that the lighter gas particles follow a 
Boltzmann velocity distribution. Starting from the kinematics of the random collisions we show that 
(f ) the average force (F) oc —x and (2) the correlation function of the fluctuating force rj = F — (F) 
is related to the strength of the average force. 
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I. INTRODUCTION 



In recent years, many experiments on single large 
molecules like colloidal spheres and bio-polymers, have 
been reported Q. Many of these experiments focus on 
the kinematics of motion of these large molecules through 
typical fluid environments, at room temperature. Most 
theoretical explanations of these experiments ^ || ne- 
glect the intrinsic correlations of the surrounding fluid 
and treat the fluid as a bath, generating friction and un- 
correlated noise on the large molecule, exactly the way a 
Brownian particle is treated. 

This approach may fail if the molecule is large or if it 
is close to a hard wall such that the displaced fluid due to 
its motion does not return to equilibrium immediately. In 
that case one may have to solve the Navier- Stokes equa- 
tion to correctly incorporate the dynamics of the fluid, 
for example while deriving the Stokes' law j^. Further, 
extended macromolecules of the polymer variety, where 
one part of the molecule disturbs the fluid which in turn 
affects the other parts of the molecule (the so called hy- 
drodynamic effect [Q) require yet another kind of analy- 
sis. It is however worthwhile and desirable to gain more 
physical insight and understanding into the simple yet 
intriguing situation of the Brownian motion, for its own 
sake and also to come up with meaningful models for 
more complicated situations. We revisit the problem of 
Brownian motion through an ideal gas in this work. 

When a particle, substantially heavier than the gas 
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particles, is injected into the gas, it is randomly kicked 
around by the gas molecules. However, due to the heavier 
nature of this particle (usually called the Brownian parti- 
cle) it takes a bunch of small kicks by the gas molecules, 
over a time interval, to move the Brownian particle ap- 
preciably. Thus there occurs a natural separation of time 
scales for the motion of this heavier particle and the gas 
molecules |^. The motion of a Brownian particle, on 
a coarse grained time scale is usually described by the 
Langevin equation j|, 



Mx + Tx = rj. 



(1.1) 
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where x is the position of the Brownian particle. The 
overall effect of collisions is modeled by introducing the 
damping coefficient T and a random force rj with zero 
mean. This is a physically motivated phenomenological 
equation where the effect of collisions is neatly separated 
into velocity dependent friction and a random force, with 
a nontrivial relationship between the two, the fluctuation 
dissipation theorcm(FDT) [||. 

In the general context of many body systems the 
Langevin equation provides a description of the slow de- 
grees of freedom, while treating the fast modes as noise 
acting on the slower ones. Apart from Brownian motion, 
thermally activated barrier crossing phenomena |^ , poly- 
mer dynamics in a solvent fluid and epitaxial growth of 
a surface profile due to random deposition are a few 
examples of such a situation and the description assumes 
that a distinct natural separation of time-scales exists be- 
tween the fast and the slow degrees of freedom. The non- 
triviality of the Langevin description lies in splitting the 
effect of the fast degrees of freedom into a systematic part 
(friction) and a random noise obeying FDT. Mostly, the 
Langevin equation is motivated from a phenomenological 
point of view. Every derivation of the Langevin equation 
requires certain assumptions; Mori and Zwanzig have 



given a formal prescription (projection operator method 
to derive the Langevin equation starting from a mi- 
croscopic picture. Given a Langevin equation, Kubo's 
method |l^ can be used to derive FDT with certain 
assumptions. Assuming the motion of the Brownian par- 
ticle through an ideal gas to be a jump Markov process, 
Gillespie has derived the Langevin equation 1 14 . A more 



physical approach to derive Langevin equation and FDT, 
though in one dimension, has also been reported [ p^ . 

In this paper we present a new approach to derive the 
Langevin equation for a Brownian particle moving in an 
ideal gas at a fixed temperature, using a combination of 
microscopic and statistical ideas. The work we present 
here differs from others mainly in the way we derive the 
FDT and that we consider two different limits of the ideal 
gas, namely, small mean free path and large mean free 
path. 

The fact that for a moving Brownian particle the rel- 
ative velocity of the gas particles approaching it from 
the front is higher than the ones approaching from be- 
hind, causes retardation leading to net damping. So the 
coUisional force, though rando m i n nature has a finite 
average (F) (= — Fi, in Eqn. (LI)). If we separate out 
this time independent average from the total force, what 
remains (i. e, F — (F)) is a time varying random force 77. 
For the Langevin description to be valid the probability 
distribution of rj must satisfy FDT. 

We explicitly show this for two different limits of the 
ideal gas, namely the one with high collision rate and the 
other with negligible collision rate. 

The material in this paper is arranged as follows: In 
section ^ we calculate the damping constant and derive 
the fluctuation dissipation theorem for a gas with a low 
collision rate. The gas molecules in this case do not suffer 
too many collisions among themselves in the time taken 
by the Brownian particle to gain a sufficient amount of 
momentum. Section [II deals with similar questions for 



an ideal gas with a high collision rate and in this case 
the gas particles can equilibrate in a short time com- 
pared to the time taken by the Brownian particle to gain 
momentum. Section IV contains a discussion of the re- 
sults and some conclusions. In appendix ^ we discuss 
the case when the Brownian particle is assumed to be a 
sphere instead of a disk and appendix ^ has some useful 
mathematical formulae 




FIG. 1: A disc shaped Brownian particle will collide with all 
the molecules with a given velocity v in the slanted cylinder 
in unit time. 



case without getting lost in the geometrical factors. The 
easiest way to arrive at the force is to treat all collisions 
in the rest frame of the Brownian particle. If a molecule 
comes in with a velocity v then from the conservation of 
energy and momentum, it follows that it is reflected with 
a velocity — kv 



M ■ 



M - 



(2.1) 



M and m being the masses of the Brownian particle and 
of the gas molecule respectively. For M ^ to, the mo- 
mentum transfered to B per collision is then 



(Sp = (1 + k)tov 



(2.2) 



The force on B is the momentum transferred to it per 
unit time, for which we need to know the number of col- 
lisions per unit time for each value of v. This number 
for a given v is the number of molecules in the (slanted) 
cylinder of edge length |v| swept in unit time as shown 
in Figure |l|, i.e. 



II. BROWNIAN MOTION IN AN IDEAL GAS 
WITH MEAN FREE PATH A ^ 00 



Consider a Brownian particle 'B' moving in a gas along 
the X direction with a velocity uo. To find the friction 
coefficient F, we have to calculate the net force on B 
due to collisions with the gas molecules. Consider the 
simplest possible geometry for B, that of a plate moving 
with a velocity perpendicular to its plane (In appendix 
the calculation is extended to a spherical particle). It is 
easier to capture the essentials of the calculation in this 
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(2.3) 



where A is the area of cross section of the plate and p 
the density of the gas. The force on B is Ncoi times the 
momentum transfer in each collision, integrated over all 
possible velocities: 



¥^-Ap{l + n) j dvPe(v)|t;,|v. 



(2.4) 



where Pc (v) is the distribution of velocities in the frame 
of the Brownian particle. The distribution in the lab 
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frame is Maxwellian i.e. 



PoM = 



27rcr 



-e 2a 



(2.5) 



where a is the mean thermal velocity component in one 
direction. 
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(2.6) 



If the velocity of a gas molecule is v in the co-moving 
frame, then it is v + uqx in the lab frame. Therefore the 
distribution in the co-moving frame is, 

Pc {Vx,Vy,Vz) = PoiVx +Uo)Po{Vy)Po[Vz) (2.7) 

The component of the force perpendicular to x adds up 
to zero whereas the x component is 



^Ap{l + K) I dvxPo{vx)\vx - uo\{vx ~ uo) (2.. 



After a few simple manipulations, we get: 

Fx = ~2Ap[l + K)e J dvxV^—j=e 2^ sinh — —j 


(2.9) 

It may be noted that the force is an odd function of uq, in 
keeping with the general rule that dissipative terms in an 
equation of motion violate time reversal invariance. For 
small velocities uq, we can Taylor-expand to the lowest 
order in uq to get 



Fx = -2Ap{l + k)uo(t- 
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(2.10) 



We call the dimensionless integral in the above equation 
X3 (Appendix ^ . The final expression for the force then 
is: 

Fx = -2Ap{l + K)uoax3 (2.11) 
from which we read off the friction coefficient 



r = 2Apil + K,)ax3 



(2.12) 



The friction coefRcient is not proportional to the linear 
dimension of the particle(as it is in the Stokes' Law) but 
rather to the area. This result carries over to the case of 
the sphere as well. The area comes because the gas we 
are considering is ideal and the motion of the Brownian 
particle in no way affects the flow of the gas. The Stokes' 



result essentially derives from the fact that the fluid (be- 
ing a strongly interacting gas) makes way for the particle 
as it moves through it |^ . 

Let us now proceed with the calc ulat ion of the noise. 
In the standard Langevin equation (1.1) the noise 77 has 
to satisfy the fluctuation dissipation theorem 



{V{t)v{t)) = 2kBTTd{t - t) 



(2.13) 



This theorem relics on the requirement that the Brown- 
ian particle is in thermal equilibrium with the gas it is 
moving through. 

Before we derive the fluctuation dissipation theorem, 
we coarse grain the time in terms of a small time unit 
e, chosen such that several collisions take place within e 
time, while the Brownian particle does not move appre- 
ciably during this interval. For a massive enough particle 
it is always possible to find such a time. It is simplest to 
calculate the noise in the rest frame of the particle , wh ere 
the discretized version of the Langevin equation (LI) is 



M 



Av 



77. 



(2.14) 



The desirable fluctuation dissipation theorem that the 
noise is expected to satisfy then reads in discrete form as 



(2.15) 



Note, that the function 5{ti ~ tj) 5ij/e, such that 
for i = j, (5(0) 1/e, which is the appropriate discrete 
form p^ . This is the relation we will now derive from 
the kinematics of the collisions between the Brownian 
particle and the gas molecules. The momentum along a;, 
transferred to B in the interval e is ery, which we call 5px- 
We have to show that: 

• The momentum transfered 5px is distributed inde- 
pendently in each time interval e. 

• The mean square value of 5px in any given time is 
given by 



{{5px) 



2kBTTe 



(2.16) 



The first of the above statements follows from the fact 
that the gas is ideal, so that each collision is independent 
and the coarse graining over time will obviously not gen- 
erate correlations. Suppose a total of N collisions take 
place in the interval e, counting all velocities. The total 
momentum transferred is |]l6| 

Spx = m{l + k){vxi + .... -I- Vxn) (2.17) 

Once again, since the velocities are independent, all cross 
terms have a zero average, leading to 



{{5pxf)^m\l + KfN{vx') 



(2.18) 



The number of particles of velocity v that collide in time 
e is 



Ar(v) = JLAePo{v)\vx\ 
m 



(2.19) 
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The total number of collisions is most easily calculated 
in Cartesian coordinates: 



N = A-^e [ dvP,(v)|w, 

TO J 



2A^aexi+0{{uo/<jf) 



(2.20) 



A calculation in polar coordinates to the leading order 
gives 



N = A^a€X3 



(2.21) 



which is the same result as Xs/xi — 2 (as shown in Ap- 
pendix^). The average {vx'^) is not the thermal average, 
but the average over the sample of particles that collide 
with B in the interval e. 



JdvN{v) 



N TO 



dvPc{v)vl\vj; 



2a' 



(2.22) 



Note that the mean square velocity in Eqn. (2.2S) is 
larger than the thermal average velocity a^. This is due 
to the fact that in the ensemble of molecules that collide 
in a given interval, there are a greater number of faster 
molecules, as they sweep a larger cylin der ( Figure |^). 
Substituting for N and (w^) in equation ( 2.18 ) we get 



{{6px) ) = m^{l + Ky2A—ecrxi2a^ 



(2.23) 



Com binin g equ atio ns ( ^.23| ) and (2.12) and using equa- 
tions (2^) and (B2) the above equation becomes 



{l + K)kBTTe 



(2.24) 



which for K = 1, i.e. for M » to goes to the desired 
fluctuation dissipation result ( 2.16| ). Also note that the 
change of momentum (or velocity) of the Brownian par- 
ticle, in a time step e, is of the order ^/e and not e, which 
might be naively expected. This is an essential feature 
of Langevin equation ||l^ . 





FIG. 2: Plot of the distribution N{v) for the ideal gas in two 
different limits. The upper plot is for an ideal gas with large 
mean free path while the lower one is for a gas with very small 
mean free path. The dotted curves in both the cases are for 
a Brownian particle at rest and the solid lines for a moving 
Brownian particle. 



III. BROWNIAN MOTION IN AN IDEAL GAS 
WITH SMALL MEAN FREE PATH 

The classical ideal gas is a familiar concept, described 
by the equation of state 

PV = nRT. (3.1) 

This equation of state is a limiting case of a more phys- 
ically reasonable equation of state, namely the real gas 
(Van der Waals gas) equation (Q) 

iP-^)iV-b)=nRT. (3.2) 



where a and b represent the intermolecular forces and 
the volume occupied by the gas molecules respectively. 
The physical meaning of working in the ideal gas limit 
corresponds to neglecting all forces between the particles 
other than the two body hard sphere collisions, and ne- 
glecting the volume occupied by the molecules themselves 
as compared to the volume of the gas. 

It turns out that there is no unique way of reaching 
such a limit and different limiting procedures lead to dif- 
ferent physical situations. One way is to make the gas 
more and more dilute in the hope that it will finally be- 
come ideal. As a consequence of dilution, the collision 
rate is reduced and finally goes to zero while the mean 
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free path approaches infinity. Though such a gas will sat- 
isfy the ideal gas equation of state, it will have infinite 
relaxation time, and will never equilibrate on its own! 
However, once the equilibrium is established it is main- 
tained in the absence of external perturbations. This is 
the kind of gas we considered in section |l| for the passage 
of a Brownian particle through it. 

The second way is more interesting, in that one can 
consistently carry out the limiting process from a real 
gas to an ideal gas while retaining a finite collision rate. 
In fact it is possible to fix the collision rate at any de- 
sired value and go through the limiting procedure. As a 
consequence we can make the relaxation time as small as 
we wish. Therefore, if disturbed from equilibrium, such 
a gas will equilibrate on its own. 

Both these cases correspond to different physical situ- 
ations and the behaviour of a Brownian particle injected 
in either of the gases is very different in each case. 

Going back to the motion of a Brownian particle in sec- 
ond type of gas we observe that the relaxation time(the 
time in which the gas goes back to equilibrium, after be- 
ing disturbed) can be made smaller than the coarse grain- 
ing time (the time for which the Brownian particle does 
not change its velocity appreciably). In fact we are al- 
lowed to assume any preassigned value for the relaxation 
time and other parameters can be adjusted accordingly. 
The picture of the previous case where we assumed that 
the gas is practically collisionless(at least within coarse 
graining time e) is no longer valid. The number of gas 
particles colliding with the Brownian particle per unit 
time at velocity v is no longer the number of particles 
in the slanted cylinder of height Vx , as the particles keep 
colliding and forgetting their velocities within the time e. 

It turns out that it is simpler to compute the colli- 
sion rate in this case as we do not have to worry about 
the slanted cylinder. The gas particles colliding with the 
Brownian particle are chosen randomly from a continu- 
ously equilibrating velocity distribution, as seen from the 
rest frame of the Brownian particle. The memory effect 
present in the slanted cylinder description of the previ- 
ous section disappears. The number of collisions per unit 
time A''(v) of the gas molecules moving with a velocity v 
with the Brownian particle in this case (Figure H) is 



m 



(3.3) 



where Vth = y ^ is the thermal velocity, K is a nor- 
malisation constant and all other parameters have the 
same meaning as in eq. (2.19) The average force can now 



be readily obtained as in the previous case 



N{v)Spdv. 



(3.4) 



Using the value of 6p from the equation (2.2) and of N{v) 
from eqn. (3.3) we have 



F ^ K{l + K)pAvthJ Pc(v)vdv 

= -uoT (3.5) 



where the friction coefficient T = k{l + K)pAvth- It is 
interesting that the friction being proportional to the ve- 
locity of the Brownian particle is an exact result here. 

We now move on to the noise calculation, the treat- 
ment being very similar to that of section The momen- 
tum transfered in each time e is obviously independent 
and we basical ly on ly have to show that (p^) = 2kBTTe. 
The equation ( ^.1§| ) is still valid for the Brownian plate 
where only the x component of momentum is important 
but the expressions for the total number of collisions N 
and (w^) are different in this case. We have for the Brow- 
nian particle at rest 



N 



N{v)dv 



K^Avth 
m 



(3.6) 



and since there are no memory effects we have the pure 
thermal velocity distribution with 



(3.7) 



Using these expressions and the equation ( ^.1§| ) we get 



{pD = m'il + nfK^Av.^e'^ 
m m 

= {l + K)ekBTT 



(3.8) 



In the limit when Brownian particle is of large mass com- 
pared to the gas molecules i.e. n — > 1 we recover the 
fluctuation dissipation theorem (2.16). The derivation 



has been simpler and more exact in this case, where the 
gas is constantly equilibrating. We have not used any 
approximations here as opposed to Section II. It is also 
clear in this case that the separation of time scales for 
the motion of the gas molecules and that of the Brow- 
nian particle is essential for the fluctuation dissipation 
theorem to hold and thus the theorem is valid only in 
the limit of k — > 1. 



IV. CONCLUDING REMARKS 

We have shown, how the familiar result of Langevin 
equation and fluctuation dissipation theorem can be de- 
rived and understood in a special case when the fluid 
through which the Brownian particle moves is an ideal 
gas. The frictional force comes about because of the fact 
that if B moves in a particular direction then as seen 
from the frame of B the velocity of gas molecules hit- 
ting it from the front is more compared to that for the 
molecules hitting it from behind. As the details of the 
calculations in section II and III show, the situation is 
a little more subtle and we have to be careful while set- 
ting up the expression for the collision rate as a function 
of velocity. Moreover, the separation of time scales for 
the motion of the Brownian particle and the gas parti- 
cles is essential for obtaining the fluctuation dissipation 
theorem. 
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As a matter of contrast we observe that the essential 
difference in the two kind of ideal gases is in the distribu- 
tion A^(v) of collision rate is at velocity v. For the ideal 
gas of the first kind, this function tends to have larger val- 
ues for larger velocities as compared to the Maxwellian, 
as the volume of the cylinder from which all the molecules 
must collide with B is directly proportional to | v| . On the 
other hand, for the ideal gas of the second type, the ab- 
sence of memory effects and the breakdown of the cylin- 
der picture prevents any such biases. 

Another interesting observation is that the friction co- 
efficient is proportional to the cross-section area for both 
the cases of the plate and the sphere unlike in the case of 
Stokes' law where it is proportional to the radius. This is 
due to the fact that we have neglected the back reaction 
due to the motion of the sphere on the velocity distribu- 
tion of the gas particles, which has been retained to be 
Maxwellian all along. On the contrary for the calcula- 
tion of the Stokes' law one takes into account the effect 
of the motion of the sphere on the fluid velocity. While 
our treatment holds for a dilute noninteracting gas, the 
Stokes' law is obeyed in fluids with higher density. 

We feel that our analysis can be carried out for other 
kind of Brownian motion as well. One could, for exam- 
ple, consider the Brownian motion of angular variables 
for molecules moving in an ideal gas . We expect that 
interesting physical insights can be gained if such an anal- 
ysis is carried out. 



APPENDIX A: CALCULATION OF FRICTION 
COEFFICIENT FOR A BROWNIAN SPHERE 
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FIG. 3: A pictorial description of collisions for a spherical 
Brownian particle. 



The velocity distribution of equation (2/7) breaks the 
0' symmetry but preserves the symmetry, therefore in 
the frame F the x and y component of the force averages 
to zero and we need to consider only the z component of 
the force. 



In this appendix, we demonstrate that the friction coef- 
ficient can be calculated for the case when the Brownian 
particle is a sphere instead of a plate. The underlying 
physics is same though the calculation is a little more 
involved. The plate was chosen only for simplicity and 
this calculation is given here to convince the reader that 
there is nothing special about the plate. We will see that 
the only difference from the plate case is a geometrical 
factor. 

Imagine the Brownian particle to be a sphere of radius 
R moving in the z direction with velocity Mo(As shown 
in Figure ^ . We consider a frame of reference F which is 
at rest with respect to the Brownian particle. Consider 
the collision of the Brownian particle with a beam of 
particles moving with velocity V(Vo, , <!>') in this frame. 
Further consider a frame F' again at rest with respect to 
the Brownian particle but with its orientation such that 
the z' axis is in the direction of V. When a narrow beam 
hits the sphere at a point {R, 9, (f) in the frame F', it 
is refiected by the sphere and the force exerted on the 
sphere, which is along z' direction is given by 

Fz'-rcst ^ fluX.(l + K).V;C0s6lz'. 

Where 

flux = pY.{-ds)^ pVo cos eR^saieded<t>{kl) 



-F* = F5,/_rrst,. cos^'. 



(A2) 



Rewriting the eqn (2.7) for the velocity distribution in 
the spherical polar coordinates as seen from the frame F 
we have 



^c(V) 



1 



[(Vo cos 0'+Uo)^ + (yo sin e')^]/2cr^ 



(V2^(7)3 



(A3) 



The total frictional force can now be readily obtained 
by integrating over all the angular variables and the ve- 
locity distribution 



JlUaVo cose' 1 2a'^ -nl 



(V27ra)3 



sin 61' cos e'dAi) 



where, 



Ci (27r)2(l+K)i?2e^ and C2 = / ^6* sin 61 cos^ 6* = 2/3 

Jo 

(A5) 

To the lowest order, the J dd' gives UoVo/3cr^ and by 
working out the / dVo we obtain the total frictional force 



F. 



DTT 



(A6) 
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Using the dimensionless integral ^5 from appendix r ^2 

the relation X5 = 4x3 and substituting the values of =2" dxx"e^ (Bl) 

Ci, C2, the final expression for the force simplifies to 



00 



^ -I {2Ap{l + n)uoax3) (A7) x« = ^ / x" e""^ = ^2^^ T ' " ^ 





where A = nR"^ is the cross section of the sphere. We Xs/xi — 2 Xs/Xs = 4 (B2) 

notice that the only difference in this expression as com- 



pared to the plate case of eqn. (2.11) is the geometrical 
factor of 8/9. 



APPENDIX B: SOME USEFUL INTEGRALS 

We list some useful integrals used in the paper. 



Tin) = / dx a;"-^ e"^ 2 / dx a;^"-! e"^ 
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